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By matching Wilson loop ratios on anisotropic lattices we measure the coefficients Ca and Cr, which are required 
for the calculation of the energy density. The results are compared to that of an indirect method of determination. 
We find similar behaviour, the differences are attributed to different discretization errors. 



1. Different ways to calculate the energy 
density 

From statistical mechanics we know the energy 
density as e = ( d(i/T) ) ■ lattice calcu- 
lations the straightforward approach to calculate 
this quantity is to introduce different lattice spac- 
ings aa, a-r to perform the derivatives with re- 
spect to 1/T = Nrttr at fixed physical volume 
{NcraaY . Therefore one introduces the lattice an- 
isotropy ^ — ^.a / ar and two gauge couplings in 
the action, 
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The energy density e contains then the derivatives 
of these gauge couplings. 
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where 

5*0 - action on a symmetric lattice, T=0 contribution, 
Scr - spatial part of the action 1 finite temperature 
Sr - temporal part of the action J simulation, 
and 
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The coefficents Ca and Cr were calculated pertu- 
batively in ||l| . 

*Talk presented by T. Scheideler, work supported by DFG- 
grant Pe 340/3-3. 



A different way to determine the energy den- 
sity was used by the Bielefeld group The free 
energy density / is obtained by integrating the 
difference of the plaquettes. 
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In large systems, where p = — /, e can be found 
from the pressure, 
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= \2NN^{c„ + Cr) [25-0 - (5, + Sr)] , 

Using the relation c„ + Cr ^ (°'§a) ^'^'^ 
/3-function from Tc/^/a measurements, we can 
calculate and Cr from (||), (||) and (||). In Fig- 
ure 1 the results for Ca and Cr are shown as solid 
lines. They deviate clearly from the pertubative 
results (broken lines). However, this is expected, 
as also Ca + Cr deviate from the pertubative re- 
sult. We also note, that the method works only 
above /3criticai(-^r), becausc both p and e become 
very small below the critical point. 

2. Direct determination of c^ and Cr by 
matching of Wilson loop ratios 

An attempt to determine the lattice anisotropy 
was made by Burgers et al.[^. Consider the an- 
isotropic couplings 
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Fig. 1. The derivatives and Cr for the Wilson 
action (a) (iV^=4,6,8) and the 2x2-action (b) (7Vr=4) 
1^ versus 13. The broken hnes show the pertubative 
values Q. The sum Co- + Cr is from the /3- function. 
The dotted line shows an estimate form the different 
TVr-values. The single measurements are from the 
matching method. 

of the action (|l|). The derivatives with respect to 
the lattice anisotropy ^ — a„jar are 
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We vary the anisotropy of the couplings 7 in the 
action S — P / ^ Sa + (3"f St and measure the an- 
isotropy of the lattice spacing ^. Wilson loops 
are suitable observables, since they depend on the 
physical size of the lattice. Instead of matching 
Wilson loops, measured in spatial and in tempo- 
ral direction, directly. 
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we match ratios of Wilson loops, in order to can- 
cel corner and selfmass contributions, 
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Since these ratios contain only Wilson loops with 
the same extension in t, the following matching 
condition holds. 
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which requires the function 7(^) or ^(7). 



Fig. 2. The ratios Ri^t{x, t) for fixed x = 1, . . . , 7 vs. 
t + x, the open symbols are connected by solid lines. 
The ratios Ri,a(x, z) correspond to the bold face sym- 
bols and are shifted in t by ^ = 1.63(4). 

We are using a 16" lattice for 7 e [0.92, 1.08], 
a 16^ X 32 lattice for 7 e [1.1, 2.0] and a 16^ x 48 
lattice for 7 = 3.0. The link integration technique 



of ref. is applied to obtain accurate expectation 
values of large Wilson loops with a statistics of 
2000 measurements. 

The value of ^ was chosen in such a way, that 
the square deviation of Ri,a{x,z — £^t) from the 
lines connecting the Ri^r{x,t) measurements was 
minimal. For each 7 we find in this way a value 
of ^. Figure 3 shows an example for f3 — 6.3. In 
all cases we found a linear behaviour of ^(7). 
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Fig. 3. The result for ^ as a function of 7 for (3 = 6.3. 
The solid line is a linear fit to ^(7). 

To obtain the derivative d£,/dj\^^-^^ we have 
performed linear fits. The resulting Ca- are shown 
in Figure 4. 

Obviously the matching procedure requires at 
least loops of an area eight. For larger Wilson 
loops the results are consistent. 

3. Summary and conclusion 

Employing a method of matching ratios of Wil- 
son loops on anisotropic lattices we have directly 
determined the derivatives and Ct of the gauge 
couplings non-pertubatively. This has been done 
for both, the standard Wilson action and the 2x2 
improved action jj]. We found significant devia- 
tions from the pertubative result. Qualitatively 
our results are in agreement with the indirect 
method used inB, though they are a little lower. 



We assume that the difference is due to ©(a") 
corrections, since they enter the two methods dif- 
ferently. 

The use of ratios of Wilson loops was crucial in 
order to eliminate unphysical self energy contri- 
butions. The numerical matching technique for 
these ratios enabled us to obtain the required ac- 
curacy for f . 
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Fig. 4. The derivative Ca for /? = 6.3 from Eq.(|7|) as a 
function of the minimal area A of the smallest Wilson 
loops included in the ratio matching. The solid line 
shows the average value for A > 8, the broken lines 
the error band. 
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